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Employing the perturbative treatment of gravitational clustering, we discuss possible effects of 
primordial non-Gaussianity on the matter power spectrum. As gravitational clustering develops, 
the coupling between different Fourier modes of density perturbations becomes important and the 
primordial non-Gaussianity which intrinsically possesses a non-trivial mode-correlation can affect 
the late-time evolution of the power spectrum. We quantitatively estimate the non-Gaussian effect 
on power spectrum from the perturbation theory. The potential impact on the cosmological param- 
eter estimation using the power spectrum are investigated based on the Fisher-matrix formalism. 
In addition, on the basis of the local biasing prescription, non-Gaussian effects on the galaxy power 
spectrum are considered, showing that the scale-dependent biasing arises from a local-type primor- 
dial non-Gaussianity. On the other hand, an equilateral-type non-Gaussianity does not induce such 
scale-dependence because of weaker mode-correlations between small and large Fourier modes. 

PACS numbers: 98.80.-k 



I. INTRODUCTION 

Recent observations of cosmic microwave background (CMB) anisotropics as well as density perturbations in the 
large-scale structure strongly support the basic predictions of inflationary scenarios, in which primordial adiabatic 
fluctuations were produced during the accelerated phase of the cosmic expansion and their statistical properties are 
approximately described by Gaussian statistics with a nearly scale-invariant power spectrum (e.g., [33l l6o|). Some 
specific inflationary models have been ruled out by observations, narrowing the constraints on the early stage of 
the universe. With precision measurements from future observations, we will detect clear signals that help us to 
discriminate between many candidate of inflationary models. 

Amongst several signals accessible in the near future, departure from Gaussianity is an important clue to probe 
the generation mechanism for primordial perturbations as well as to discriminate between inflation models. While 
the simplest slow-roll inflation with single scalar field predicts a small departure from Gaussianity (e.g., pi [Til [4l| ) . 
models with a non-trivial kinetic term or late-time inflaton decay called curvaton scenario [37l . |48| can produce a large 
non-Gaussianity (e.g., [j| [H, HH, HI] ) . There are also viable models motivated by string theory that pro duce large non- 
Gaussianity such as Dirac-Born-Infeld (DBI) inflation 0, @, 0, GJ, HE HE IH| and ekyprotic scenarioMjp, [H, |M [H . 
Although tentative detections of primordial non-Gaussianity have been reported very recently f |47l. l69j. but see also 
[HI, HH) and the result is broadly consistent with the standard prediction of slow-roll inflation, these detections are 
at relatively low statistical significance and more precise measurements are necessary for a definite detection of the 
non-Gaussian signals. Planned CMB surveys, namely the Planck mission (52|, will have much greater sensitivity 
and better detect non-Gaussian signals (e.g., |6l|). Further, large-scale structure data from future spectroscopic 
surveys will uncover the mass density fluctuations up to Giga parsec scales, which preserve the statistical properties 
of primordial fluctuations (e.g., [57L [5^]). Since these two measurements can probe different scales of primordial 
fluctuations, they are, in principle, complementary to each other. 

In the present paper, we study the signature of primordial non-Gaussianity imprinted on the large-scale structure, 
especially focusing on the matter power spectrum. Naively, we expect that the signature of primordial non-Gaussianity 
basically appears in the statistical properties of higher-order quantities such as the bispectrum and trispectrum, and 
the power spectrum as a second-order statistic remains unchanged even if large non-Gaussianity has been produced. 
However, this is true only when the fluctuations of the mass density field are tiny and well within the linear regime. 
As the gravitational clustering develops, the coupling between different Fourier modes becomes important and the 
scale-dependent non-linear growth appears due to the mode-coupling effect. In the weakly non-linear regime, the 
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strength of this mode-coupling sensitively depends on the initial condition. Since the non-Gaussian density field 
intrinsically possesses non-trivial mode-correlations, the evolution of the power spectrum may be altered by the 
primordial non-Gaussianity, at least in the weakly non-linear regime. 

Here, employing the perturbative calculations, we study the effect of primordial non-Gaussianity on the matter 
power spectrum. The effect of non-Gaussianity on the matter power spectrum has been previously investigated by 
Ref. [H3| using perturbation theory, showing that the magnitude of this effect is roughly at the few percent level in 
the weakly non-linear regime. In this paper, adopting the two representative models of primordial non-Gaussianity 
described in Sec. [TTJ we quantitatively estimate the non-Gaussian effects on the power spectrum and elucidate their 
shape dependence (Sec. IIIip . In particular, we carefully examine the influence of primordial non-Gaussianity on 
cosmological parameter estimation, paying particular attention to the primordial spectral indices and the characteristic 
scale of the baryon acoustic oscillations as a cosmic standard ruler fSec. IIV|) . Also, we discuss the non-Gaussian effect 
on the galaxy biasing (Sec. JVJ . This issue has been recently studied by several authors in the context of the halo 
biasing prescription [2. [T9L [42L l4rl |62| . and local biasing prescription Here, we adopt the local biasing scheme 
and consider the scale-dependent galaxy biasing arising from the primordial non-Gaussianity. 

Throughout the paper, we assume a flat Lambda cold dark matter (CDM) model and the fiducial model parameters 
are chosen based on the five-year WMAP results (WMAP+BAO+SN ML, see Ref.[H|): Q h h 2 = 0.02263, tt c h 2 = 
0.1136, fl K = (fi A = 0.724), n s = 0.961, r = 0.080, A n {k = 0.002Mpc _1 ) = 2.42 x 10~ 9 , h = 0.703. 



II. PRIMORDIAL NON-GAUSSIANITY 



with the function Mq being 



According to the commonly accepted mechanisms to generate the primordial fluctuations, non-Gaussianity would 
be imprinted on the primordial curvature perturbation produced during inflation. The curvature perturbation can also 
be generated after inflation by a late-time decay of light fields, and this would produce a large non-Gaussianity in the 
curvature perturbation. Further, as alternative scenarios to inflation such as ekpyrotic models, nearly scale-invariant 
curvature perturbation with large non-Gaussianity could be generated during the collapsing phase of the Universe. 

Let us denote the primordial curvature perturbation on the uniform density hypersurface and the Bardeen's cur- 
vature potential at super-horizon scales by Q p and &h, p , respectively. The primordial density perturbation defined at 
an early time of the matter dominated epoch, So, is related to these variables through (e.g., [25L I39I IBH l59l|) 

6 (k) = M c (k)( p (k) = ^M c (k)<Z> H , p (k), (2.1) 

M ^ = lP%> (2 - 2) 

where the function T(k) is the transfer function of matter fluctuations normalized to unity at k — » 0. Note that 
the curvature potential given in the matter dominated epoch, is related to the curvature perturbation C(fc) = 

(5/3) $ Hj ,(fc). 

Using the relation (|2.1|) , n-point statistics of the primordial density field are all expressed in terms of respective 
correlator of the curvature potential or the curvature perturbation. For the power spectra, we have 

P (fc) = M 2 (k)P c (k), (2.3) 

with quantities Pq and P^ defined by 

(S (k)So(k')) = (2ir) 3 5 B (k + k') P (k), (( p (k)( p (k')) = (2vr) 3 &(fc + k') P ( (k), (2.4) 
The bispectrum as the first non-trivial quantity arising from the non-Gaussianity becomes 

Bo(k 1 ,k 2 ,k 3 )=M c (k 1 )M c (k2)M ( (k 3 ) B ( (k u k2,k 3 ), (2.5) 

where we define 

(<5o(fci)<5o(fc2)5 (fc3)) = (2^) 3 S D (ki + k 2 + k 3 )B (k u k 2 ,k 3 ), 

|C P (fci)C P (fc 2 )C P (fc 3 )) - (2^) 3 ( 5 D (fe 1 +fc 2 + fc3)B c (fc 1 ,fe 2 ,fc 3 ). (2.6) 

Details of the functional form of Bq depend on the generation mechanism of primordial non-Gaussianity. In the 
following, we will consider two representative models of non-Gaussianity and give the explicit expressions for the 
primordial bispectrum Bq. 
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A. Local model 

In the local model of primordial non-Gaussianity, the primordial fluctuation <&H,p is characterized by the Taylor 
expansion of the Gaussian field. Denoting the Gaussian field by <p, it is conventionally characterized as: 

$ H , p (x) = <p(x) + / NL {v 2 {x) - (<P 2 )} +■■■ ■ (2.7) 
The Fourier counterpart of the above equation is given by 

* ff)P (fc) - <p(k) + /nl J ^3 MqMk -q)- (<p(q)<p(k - q))} + ■■■. (2.8) 

Alternatively, we can expand ( p as: 

C P (fe) - Cc(fe) + ^/nl J {C G (q)CG(fc - q) ~ (CG(q)CG(fe - q))} + • • • , (2.9) 

where (q is Gaussian field and we have the relation £q = (5/3)<£>. For the cases of interest here, deviation from 
Gaussianity is generally small and perturbative evaluation of the primordial bispectrum is valid. The leading-order 
result becomes 

B c {k u k 2 ,k 3 ) ~ ~/ NL [Pdk^P^h) + P c (k 2 )P c (k 3 ) + P^P^k^. (2.10) 

In the slow-roll inflation models, /nl is suppressed by the slow-roll parameters and thus it is expected to be very 
small. However, in curvaton models or ekpyrotic models, it is possible to generate a large non-Gaussianity of the local 
model with / NL - C>(10 2 ). 

B. Equilateral model 

As another type of bispectrum configuration, we consider an equilateral model. The DBI inflation is a typical 
example where the bispectrum has a maximum amplitude for equilateral configurations. This is in contrast with the 
local model where the bispectrum has the largest amplitude for squeezed configurations [?|. The bispectrum in the 
equilateral model is approximated by [13] 1 

B c (ki,k 2 ,k 3 ) - f f^ L (j^j [-Pdki)Pdk2)-Pdk2)Pc(k 3 ) - Pdh)Pdki) 

-2 {Pdk 1 )Pdk 2 )Pdh)} 2/3 + | (Pdki){Pdk2)} 2 {Pdk 3 )} 3 ) 1/3 + (5 perm.)} ] , (2.11) 

with k — (k\ + k 2 + k-i)j?>. Here, we introduce the slope index k in order to characterize the scale-dependence of 
the non-Gaussianity, which generically appears in the DBI- type inflation models, we set the characteristic scale to 
kcMB — 0.04 Mpc -1 . Note that the bispectrum in equation (|2. 1 1|) is normalized in such a way that, for equilateral 
configurations (fci — k 2 = k 3 = k), ii coincides with the local form given in equation (|2.10p . 

III. NON-LINEAR EVOLUTION OF MATTER POWER SPECTRUM FROM PERTURBATION 

THEORY 

A. Perturbation Theory 

Even if the primordial fluctuations are well described by linear theory, the non-linearity of the gravitational dynamics 
eventually dominates and we must correctly take into account the non-linear growth of matter fluctuations, which 



1 The expression given here is slightly generalized in a sense that we do not necessarily assume a power-law form of the power spectrum, 
P c (fc). 
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significantly modifies the matter power spectrum. For the scales of interest here, especially the accessible scales of 
future galaxy redshift surveys, the non-linear evolution is rather moderate and perturbative treatment is valid. The 
mass density field is described by 

8{k;z) = 5^(k;z) + 6 {2) {k;z) + 5^{k;z) + --- . (3.1) 

In the above, the function <5^^(fc; z) represents the linear fluctuation and it is given by S^(k; z) = D(z)6 (k), where 
D(z) is the linear growth rate. Neglecting the decaying mode of linear perturbation, the solutions for higher-order 
quantities are formally expressed as [12j | 

(2 ^ 3 (n-i) n S ^ k - fcl fe «) ■ " ' k rMki) ■ • ■ 6 (k n ), (3.2) 

with Fsym being the symmetrized kernel for the n-th order solution. Then, the power spectrum of the density field, 
defined as 

(S(k; z)6(k'; z)) = (2tt) 3 (5d(A; + k') P(k; z), (3.3) 
can be calculated by substituting Eq. (|3.2p into the above expression. The resultant expression is summarized as 

(3.4) 



dW(k;z) = [D(z)Y 



P(k; z) = D 2 {z) P (fc) + P (12 \k; z) + P {22 \k; z) + p( 13 >(fc; z) 
where the terms P^ 12 \ p( 22 ) and p( 13 ) are the so-called one-loop power spectra given by 



P^Hk-z) = 2D 3 {z) J A^F^ m (q,k-q)B (-k,q,k-q), 
pW(k;z) = 2D\z) J -0^{F( 2 l( q ,k-q)} 2 P o {q)P o (\k-q\) 

F syl(Pi k ~ P)F£L{l, k - <?) T o(p, k-p,q,-k- q), 



+D\z) 
P {13) (k;z) = 6P 4 (z) 
+2D 4 {z) 



d 3 q 
(2^3 
d 3 pd 3 q 

(2tt) 6 ~ sym " 

^Fgl(k, q> -q)Po(k)P (q) 



d 3 pd 3 q 



F ^L(P, Q,k-p-q) T (-k,p, q,k p q). 



(3.5) 



(3.6) 



(3.7) 



(2tt) 6 sym ' 

These are the most general expressions for the one-loop power spectra in the presence of primordial non-Gaussianity. 
The quantity P^ 12 ) is the first non-trivial correction associated with the primordial density bispectrum B${ki, &2, £3). 
While the quantities p( 22 ) and p( 13 ' represent the first leading-order corrections in the case of Gaussian initial 
conditions, additional contributions originating from the primordial density trispectrum, To(fci, k$, k%, £4), also arise. 
Thus, in general, we need the explicit functional form of the primordial trispectrum as well as primordial bispectrum 
in order to compute the one-loop power spectra. However, in most general cases, the primordial trispectrum is of 
the order of M^Pq, which is negligibly smaller than the bispectrum, Bq ~ M 3 Pq. Since the terms including Pg 
may be regarded as the two-loop order, we can safely neglect the trispectrum contribution. With this treatment, the 
terms p( 22 ) and p( 13 ) reduce to nothing but the standard one- loop spectra, and they are explicitly given by (e.g., 

[Iliiliiliil) 



P^(k;z) 
P (13) (fc;z) 



D\z) 
D\z) 



k 3 

(27T) 2 
k 3 

1 

252 x 2 



dxx 2 Pa(kx) / dfiP$(k\/ 



1 + x 2 



2fi x 



3x + 7fi - I0fi 2 x 



-1 2 



7x(l + x 2 - 2fix) 



P (fc) / dxx 2 P (kx) 



12 



158 + lOOx 2 - 42a; 4 + {x 2 - \f(lx 2 + 2) In 



x + 1 



x - 1 



(3.8) 



(3.9) 



Also, the non-trivial contribution from the primordial bispectrum is expressed in terms of the bispectrum of the 
curvature perturbation as 



P {12 \k;z) = D 3 {z) 



k 3 
'2nf 



dxx x 



+1 rf 3x + 7fi- lOfi 2 x 
1 7x{\ + x 2 — 2fi x) 



x M c (k)M c (kx)M c (ky / l + x 2 - 2 fix) B c (k, kx, ky/l + x 2 - 2 fix). 



(3.10) 
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FIG. 1: Ratio of the power spectrum to the smoothed reference spectrum for the local model (left) and the equilateral model 
with ac — (right), for z — 1. The smooth spectra are obtained from the no- wiggle approximation of the linear transfer function 
according to Ref. [20]. The lines from top to bottom respectively indicate /nl = 300 (solid), /nl = 100 (long-dashed), /nl = 
(dotted), /nl = —100 (short-dashed) and /nl = —300 (dot-dashed). 



B. Results 



Based on the expressions given above, we now calculate the matter power spectrum, focusing on the scales rel- 
evant for future galaxy surveys. Figure [T] plots the ratios of power spectra to the smooth reference spectrum, 
P(fc)/P no _ w i g gi e (Ae), given at z = 1. The reference spectrum P n o-wiggie(^) is the linear power spectra computed 
from the no- wiggle approximation of the transfer function given by Ref. [2fJ. The left panel shows the results for 
the local model, and the right panel plots the results for the equilateral model with index k — 0. The overall be- 
havior at higher k is basically dominated by the corrections coming from the standard 1-loop power spectra, p( 22 ) 
and p( 13 \ but small variation of the power spectrum amplitudes manifests at higher k. This is the contribution of 
primordial non-Gaussianity, P^ 12 \ originating from the primordial bispectrum. Depending on the sign of /nl, the 
primordial non-Gaussianity enhances or suppresses the non-linear growth of the power spectrum and the deviation 
from the Gaussian result (/nl = 0) becomes more significant as |/nl| increases. Although the effect of primordial 
non-Gaussianity seems rather mild within the current constrained values of —9 < /jJjL ~ m an< ^ — 151 ^ /nl ~ 253 
(95%GL.)[33|, we do not immediately exclude the possibility of significant impact on the cosmological parameter 
estimation. We will discuss this issue in the next section. 

Next we will focus on the shape of the power spectrum. In Figure [5J the ratios of the power spectra to those in 
the Gaussian case (/nl = 0), i.e., P(fc; /nl ^ 0)/P(k; /nl = 0), are plotted for various redshifts and values of /nl- 
Note that in computing the ratio, both the numerator and denominator are calculated from perturbation theory. In 
Figures [1] and [21 the vertical arrows labeled by and fc 3 % indicate the maximum wavenumber below which the 
precision level of perturbation theory prediction is expected to be better than 1% and 3%, respectively. According to 
recent numerical experiments, the maximum wavenumber is empirically determined through [50J 

k 2 I" ^D 2 (z)P (q)=C, (3.11) 
Jo 6vr 2 

where C = 0.18 (0.3) for k\% (k 3 %). Solving the above equation with respect to k, we obtain k\% and k 3 % as 
function of redshift. Note that the maximum wave numbers given above have been derived by comparison between 
N-body simulations and theoretical predictions, and the resultant convergence range is narrower than those previously 
proposed (e.g., [3(| [E§|). In this sense, Eq. ([3.11| may be regarded as a conservative criterion. 

Bearing the limitation of perturbation theory in mind, we see that the deviation of the power spectrum from the 
Gaussian case also depends on the redshift and scales as well as the model of primordial non-Gaussianity. Monotonic 
change of the power spectrum amplitude is broadly consistent with recent N-body simulations in the case of the local 
model [HI, although the fractional change is typically < 2 — 3%, as expected from Figured] Note, however, that 
the level of this changes roughly corresponds to the sensitivity achievable with future redshift surveys. As references, 
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FIG. 2: Ratios of power spectrum, P(k; /nl 7^ 0)/P(k; /nl = 0), for 2 = 3 (left), 2 (middle) and 1 (right). From top to bottom 
panels, the results for the local, equilateral with k = 0, 0.3, and —0.3 are shown, respectively. In each panel, we plot the 
cases with non-Gaussian parameter /nl = +300 (solid), +100 (long-dashed), —100 (short-dashed) and —300 (dot-dashed). The 
vertical arrows labeled by ki% and k 3 % indicate the maximum wave number below which the perturbation theory predictions 
are reliable with a precision of 1% and 3% level, respectively, according to the criteria (|3.11[) [5C$ . As references, the error bars 
limited by the cosmic variance are plotted in the right panel, assuming the survey volume of V B = 4 (ft _1 Gpc) 3 (thin solid) and 
V s = 10 2 {h^Gpcf (thick solid). 



we plot the error bars in each right panel, which represent the expected I — a errors, AP(k), limited by the cosmic- 
variance, i.e., AP(k)/P(k) = (2/Nk) 1 ^ 2 with Nk being the number of Fourier modes within a given bin at k. Here, 
we specifically consider the two representative cases of V s = 4(/i _1 Gpc) 3 (thin) and V s = 10 2 (ft. _1 Gpc) 3 (thick), 
corresponding to the future surveys dedicated for the measurement of baryon acoustic oscillations (BAOs) from the 
ground and space, respectively. Naively comparing the error bars to the amplitudes of the non-Gaussian effects, 
ground-based BAO surveys with a typical volume of V ~few ft,~ 3 Gpc 3 will find it difficult to detect the signature of 
primordial non-Gaussianity only from the power spectrum. However, idealistic surveys with huge volumes seem to 
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show the potential for a definite detection of non-Gaussian effects even with the currently constrained values of /nl- 



IV. IMPACT ON COSMOLOGICAL PARAMETER ESTIMATION 



A. Fisher-matrix analysis 



Apart from primordial non-Gaussianity, there arc several parameters that affect the shape of the power spectrum. 
Among these, the primordial spectral index n s and running index a = dn s /dlnk monotonically change the power 
spectrum, which resemble the effect of primordial non-Gaussianity at some wavenumbers. A natural question arises 
whether the primordial non-Gaussianity is indeed detectable or not considering the degeneracies and how the non- 
negligible effect of non-Gaussianity adversely affects the estimations of primordial spectral index and the running of 
the index. On the other hand, for the scales accessible to the future surveys, a measurement of the characteristic scale 
of BAOs is an important clue to probe the late-time acceleration of the universe e.g., [271. 143l.l60j|). and a percent-level 
determination of the acoustic scale is required for the determination of the dark energy equation of state 0, [5l[ . In 
this respect, even the small effect of primordial non-Gaussianity may affect the determination of the BAO scale. 

Here, we apply the Fisher-matrix method to address these issues and explore the potential impact on cosmological 
parameter estimation as well as the determination of the BAO scale. For this purpose, we consider a rather simplified 
setup, namely that our observable is the real-space power spectrum free from the redshift-space distortion. Under the 
assumption of linear galaxy biasing, the observed power spectrum may be written as (e.g., [30l|): 

ftb.(M = (j^\ 3 biP masB z) , (4.1) 

\^V,tmc(z) J \ V V (Z) ) 

where the power spectrum -P mass is computed from perturbation theory. Here, b\ is the linear biasing parameter, and 
Dy(z) represents the cosmological distance averaged over three-dimensional space, Dy{z) = [D 2 A (z)/H (z)] 1 / 3 . The 
subscript "true" denotes the quantities estimated from the true cosmological model. 
Following Refs.[66|, [6?|], the Fisher matrix for the galaxy power spectrum is given by 



F _ V. r max ^^ d\nP(k-z)d\nP(k-z) f n gal P(fc;z) ( 

tij ~ (27T)* A min d6i ae j { n gal p{k- z) + 1 f ' ' 4 " ' 



where 6i represents one from a set of parameters. The quantity V s is the survey volume and n ga i is the number 
density of galaxies. The range of integration [fc m ; n , k max ] should be determined through the survey properties, and in 

1 /3 

particular the minimum wave number is limited to 2ir/V s . Note that the choice of fc m ; n may be crucial for the power 
spectrum measurement on large-scales k < O.Ol/iMpc -1 , where the scale-dependent effect of galaxy biasing becomes 
prominent through the non-trivial mode-coupling from primordial non-Gaussianity 0, EH, El |46| . |62| . Here, we fix 
kmin = 0.01/iMpc" 1 and allow fc max to vary. The effect of scale-dependent biasing will be discussed in next section. 

In the Fisher-matrix analysis presented below, we include the five parameters given by Oi = 
{n Sl a, Dy I '-Dv,truei /nl, b\). The fiducial values of these parameters are set to n s — 0.961, a = 0, IV/-Dy,true = F 
and /nl = 0. We then consider a space-based BAO survey as an idealistically gigantic galaxy survey along the lines 
of SPACE and ADEPT ( @, H3, Hi] , see also Ref. [EH). We adopt the survey parameters of a space-based BAO exper- 
iment as follows: z = 1.5, V s = 100 (h^Gpc) 3 , h = 3.25, and n ga i = 3.25 x lO^/i^Mpc)" 3 (e.g., H). Below, we 
mainly focus on the local model of primordial non-Gaussianity. The results for the equilateral model are qualitatively 
the same as in the local model and are briefly summarized in Table |U 



B. Results 



Figure [3] shows the marginalized 1-er (68%C.L.) errors on the spectral index n s (left), the running of the index a 
(middle), and the distance scale Dy / Dy,true (right) as function of cutoff scale, fc max . Here, the solid lines represent 
the results assuming no prior information to these parameters (prior 1, solid). The short- and long-dashed lines 
correspond to the results taking account of the priors on n s and a, which are expected to come from the upcoming 
Planck CMB experiment : An s = 0.01 (prior 2, short-dashed); An s = 0.01 and Aa = 0.002 (prior 3, long-dashed). 
We assume Gaussian priors for these cases. 

Figure [3] indicates that the CMB priors play an important role for the determination of the primordial indices and 
distance scales. This has been repeatedly pointed out in the literature, but we emphasize that it is indeed crucial 
when we include the non-Gaussian parameter /nl in the analysis of cosmological parameter estimation. To clarify 
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FIG. 3: Predicted 1-a (68%C.L.) errors on the spectral index n s (left), running of the index a (middle), and distance scale 
Dy I Dv,true (right) as function of maximum wavenumber fc max , assuming the survey parameters of z = 1.5, V s = 100/i _3 Gpc 3 , 
bi = 3.25, and n ga \ = 10 _4 ft 3 Gpc -3 , as an illustrative example of space-based BAO missions. Here, we specifically treat the 
local model of primordial non-Gaussianity. The solid (prior 1), short-dashed (prior 2), and long-dashed (prior 3) lines represent 
the results under the different priors (see text for details). Thick lines show the one-dimensional errors marginalized over the 
four parameters (i.e., n 3 , a, Dv / Dv,tiuc, /nl), and thin lines represents the error excluding the non-Gaussian parameter, /nl. 




FIG. 4: Two-dimensional joint 68%C.L. constraints on /nl and n s (left), a (middle), and ZV/-EV,truc (right), assuming the 
local model of primordial non-Gaussianity. In deriving the constraints, we fix the maximum wavenumber to fc max = &i% — 
0.155/iMpc" 1 and adopt the survey parameters of a space-based BAO mission. The different lines indicate the results from 
imposing different prior information on n s and a: no priors (solid): An s — 0.01 (short-dashed): An s = 0.01, Aa = 0.002 
(long-dashed). 



the influences, we exclude /nl from the Fisher-matrix analysis and evaluate the marginalized errors again. The 
results are plotted in thin lines. Comparing thick with thin lines, we see that the inclusion of the extra parameter 
/nl significantly worsen the constraints on the primordial index and its running. With strong priors on n s and 
a (prior ,?), tighter constraints comparable to those ignoring primordial non-Gaussianity would be obtained. By 
contrast, the determination of distance scale Dy / Dv.tmc is insensitive to primordial non-Gaussianity if we choose 
the cutoff wavenumber fc max around fci%, labeled by the vertical arrows. This is partly because the distance scale 
is mainly determined through the BAOs at linear scales, where the effect of primordial non-Gaussianity becomes 
negligible at z a 1.5. It is interesting to note that all the curves depicted in Figure [3] become almost constant 
at k > kx% ~ 0.155/iMpc -1 . Since this scale is slightly larger than the one inferred from the mean separation of 
galaxies (n^ ~ 0.07ft.Mpc -1 ), the constancy of the predicted errors implies that the shot-noise contribution becomes 
significant at k > k\%. 

In Figure [H we next plot the two-dimensional constraints on the non-Gaussian parameter /nl and one of the free 
parameters n s , a and Dy / Dv.tmc, marginalizing over the other remaining parameters. In this plot, the cutoff scale 
is fixed to fc max = fci%, and the 1-a (68%C.L.) contours are shown for the three different priors. As expected from 
Figure [31 the primordial index n s and its running a without CMB priors (prior 1) exhibit a strong degeneracy with 
the /nl parameter, while the degeneracy between /nl and distance scale is very weak. The uncertainty in the /nl 
parameter is therefore very large in this case. The errors can be reduced significantly when the strong constraints 
on n s and a are obtained (prior 3). However, the uncertainty in /nl is still A/nl ~ 500, which is quite a bit larger 



9 



than other observational techniques using the bispectrum or cluster abundance [39j, l57|, [59J ■ The uncertainty may be 
reduced if we increase the cutoff scale k max , although it is generally difficult to accurately predict the power spectrum 
at these scales, including the higher-order corrections. Hence, it is difficult to detect the non-Gaussian signals even 
from an idealistic large- volume survey. 

Table U summarizes the results of the one-dimensional marginalized errors for various cases with different fiducial 
values, including the case of the equilateral model. Here, in addition to the priors on n s and a (prior 3), we further 
put a Gaussian prior on k with An = 0.5, the reason for which comes from the fact that k is the slow-roll parameter 
and should be restricted to k < 1. Because of the additional parameter, the resultant errors on /nl significantly 
increase for the equilateral model, especially with n > 0. By contrast, the errors onn s , a and distance scale are 
hardly affected. 

Finally, as a representative example for a more near-term project, we consider the case of a ground-based BAO 
survey. The survey parameters we adopt are z = 1, V s = 4(/i _1 Gpc) 3 , bi = 2, and n ga \ = 10~ 3 (h~ 1 Mpc)~ 3 (@, HI], 
see also [11]). Figure [5] shows the expected one-dimensional error on -Dy/£V,truo (left) and the two-dimensional errors 
on /nl and Dv/Dv,uuc (right). Compared to the large volume of a space-based BAO experiment, the statistical 
errors of the measurement of the power spectrum are somewhat larger and thereby the precision of the distance-scale 
estimation is rather degraded. Nevertheless, the results are hardly affected by the primordial non-Gaussianity. 

In conclusion, the presence of primordial non-Gaussianity potentially affects the estimation of the primordial index 
and the running, and the CMB-based priors on these parameters would be crucial, but the distance-scale measurement 
is relatively insensitive to it for the appropriate choice of fc max . In the end, the primordial non-Gaussianity is hardly 
detectable from galaxy redshift surveys even with idealistically large survey- volume. 

V. ON THE LOCAL BIASING PRESCRIPTION FOR GALAXY POWER SPECTRUM 

So far, we have mainly dealt with the matter power spectrum and the galaxy power spectrum has been only 
considered under the simplified assumption of linear galaxy biasing. In reality, the statistical relation between galaxy 
and mass is generally complicated due to the non-linear nature of galaxy formation processes, and the linear biasing 
prescription would not hold even at large scales accessible to the future galaxy surveys. 

Although the extension of the linear biasing prescription to include non-linear effects is rather non-trivial, a straight- 
forward and frequently used prescription for galaxy biasing, coupled to perturbation theory, is the local biasing [2l|, 
in which the galaxy density field S ga \ at a given position x is described as the local function of mass density field at 
the same position, i.e., 5 ga i(x) = f[S(x)]. On large scales of our interest, this can be expressed as the Taylor series 
expansion: 

6 gal (x\R) = f[6(x\R)]=h6(x\R) + ^{[6(x\R)]*-([6(x\R)]*)} 

+ |{[o>|i?)] 3 -([ >| J R)] 3 )} + •■■, (5.1) 

where the quantities S sa \(x\R) and S(x\R) are the galaxy and mass density fields smoothed over the radius R centered 
at the position x. In the above expression, the coefficients 62 and 63 describe the non-linearity of galaxy biasing, which 



TABLE I: Marginalized one-dimensional errors (68%C.L.) on n 3 , a, Dy / Dv.truc and non-Gaussian parameters /nl and n from 
Fisher-matrix analysis adopting the survey parameters of a space-based BAO experiment 



Fiducial model 








An a 


Aq 


A(LV/LV,truc) 


A/nl 




Local (/nl = 0) 








0.0075 


0.0019 


0.0026 


335 




Local (/nl = +100) 






0.0075 


0.0019 


0.0026 


338 




Equilateral (/nl 


= +250, 


K 


= -0.3) 


0.0072 


0.0019 


0.0025 


596 


0.48 (0.37) 


Equilateral (/nl 


= +250, 


K 


= 0) 


0.0086 


0.0020 


0.0025 


1306 


0.48 (0.45) 


Equilateral (/nl 


= +250, 


K 


= +0.3) 


0.0082 


0.0018 


0.0025 


2774 


0.49 (0.49) 


Equilateral (/nl 


= -150, 


K 


= -0.3) 


0.0071 


0.0019 


0.0025 


498 


0.49 (0.43) 


Equilateral (/nl 


= -150, 


K 


= 0) 


0.0085 


0.0020 


0.0025 


1230 


0.49 (0.48) 


Equilateral (/nl 


= -150, 


K 


= +0.3) 


0.0082 


0.0018 


0.0025 


2753 


0.50 (0.50) 



parentheses represent the marginalized l-cr error when further adopting the Gaussian prior of A/nl = 75 
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FIG. 5: Expected 1-a (68%C.L.) error on the distance scale as a function of maximum wavenumber (left) and two-dimensional 
joint 68% C.L. constraints on /nl of the local model and Dv / -CV.truc fixing the maximum wavenumber to fc max = ki% ~ 
0.131/iMpc -1 (right). Here, we specifically adopt the survey parameters of z = 1, V s = 4/i _3 Gpc 3 , bi = 2.0, and n ga i = 10 -3 
/i 3 Mpc -3 , as a representative example of ground-based BAO surveys. The meaning of the lines are the same as in Figures [3] 
andgl 



are incorporated into the galaxy power spectrum P ga i(fc), through the perturbative calculation of the matter power 
spectrum ([HI IH, [f53|], see also [65[ for different parametrization schemes). Although the relation (|5.ip is just a 
phcnomenological prescription, it has been recently applied to the characterization of the BAOs, and the model of 
galaxy power spectra has been tested against numerical simulations of halo/galaxy clustering in the case of Gaussian 
initial condition [3l|, [63| . 

When we consider primordial non-Gaussianity, several non-trivial corrections to the galaxy power spectrum -P ga i(fc) 
appears, which can further alter the shape of the power spectrum on large scales. This has been first pointed out 
by McDonald 46,]. Here, we rephrase his finding and advocate its observational importance, together with potential 
problems. 

Let us focus on the scales larger than the characteristic scale of BAOs, above which the higher-order perturbations 
can be safely neglected. Collecting the relevant terms in the expansion (|5.ip . the power spectrum of the galaxy density 
fields smoothed over the radius R may be expressed as 

P gal (fc; z\R) = b\ { D 2 (z)P Q (k) + P( 12 )(fc; z) } W 2 {kR) + bxb 2 F{k; z, / NL |P) + ■■■ , (5.2) 

up to the third-order in mass density field, i.e., O(Sq). Here, W(x) is the filter function defined in Fourier space. 
In equation (|5.2p , there are two types of non-Gaussian contributions. One is the term P( 12 \k;z) coming from the 
matter power spectrum, and the other is the function F(k] z, /nl) arising from the non-linear mapping of the galaxy 
biasing. The explicit expression of the function F is 

F(k;zJ NL \R) = D 3 (z) J jLjL M c (k)M ( (q)M ( (\k - q\) 

x W{kR)W(qR)W(\k - q\R) B ( (-k, q,k q). (5.3) 

The functional form of F is clearly different from the perturbative corrections of the matter power spectrum, p( 12 ) 
fEq. (|3.5p or (|3.8p ). Hence, the presence of this term potentially leads to the scale-dependent galaxy biasing, relative 
to the matter power spectrum. 

Before continuing, we note that the quantity P ga i(k; z\R) manifestly depends on the smoothing radius R. In the 
usual sense, P ga ,i(k; z\R) would not be directly related to the power spectrum derived from the unfiltered observations. 
It has been put forward in Ref. [63| that the unfiltered galaxy power spectrum P ga i(k; z) can be recovered through 
the simple operation: 



(5.4) 
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While the above definition still includes the smoothing radius, we regard equation (|5.4p as a direct observable and 
evaluate it with the appropriate smoothing radius to see what happens. Below, we will separately treat the galaxy 
power spectra for the local and equilateral models of primordial non-Gaussianity. 




10" 2 10-' 10- 2 10-' io-2 io- 1 
k [h Mpc- 1 ] k [h Mpc- 1 ] k [h Mpc" 1 ] 



FIG. 6: Ratio of biasing factor, b(k; z)/bc(k; z) given at z = 1, in the case of the local model. Gaussian smoothing is adopted 
in order to compute the biased power spectra. R — l/i _1 Mpc (left); R = 5/i _1 Mpc (middle); R — 10/i _1 Mpc (right) 
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FIG. 7: Ratio of biasing factor, b(k; z) /ba(k; z), given at z = 1 in the case of the equilateral model. In this plot, we do not 
consider Gaussian smoothing, k = (left); k — 0.3 (middle); k, — —0.3 (right). 



A. Local model 



Figures [S] plots the ratio of the biasing factor, b(k; z)/b G (k; z), plotted against the wavenumber, in the case of local 
model. Here, the biasing factor b(k; z) is defined by 
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The function b G (k; z) is similarly defined as above, but with /nl = 0. For illustrative purposes, we specifically set the 
biasing coefficients to b\ — 2 and b 2 — 0.5, and the results for z — 1 are shown for different values of non-Gaussian 
parameter /nl- Three different panels show the dependence of the smoothing radius, adopting Gaussian smoothing, 
W{x) = exp(-x 2 /2), with radii R = l/i _1 Mpc (left), 5/i _1 Mpc (middle), and lO/i^Mpc (right). 

Clearly, the ratios of the biasing factor exhibit a strong scale-dependence, especially at k < fJ.01ft.Mpc~ 1 . Note that 
the function F logarithmically diverges unless we introduce a large smoothing radius R. Accordingly, the strength of 
the scale-dependence sensitively depends on R as well as /nl- The origin of this scale-dependence can be deduced 
from a simple manipulation as follows. Substituting the bispectrum of the local model (|2.10p into the expression (|5.2[) . 
we take the limit k — > 0. The correction P^ 12 \k) becomes negligibly small, and the term including the function F is 
the only dominant contribution to the galaxy power spectrum. We have 



F(k; Z , /nl I R) * f /nl D»(*) Po(fe) W{k *£{ R) : 



a 2 (R) = J 



^L W 2 ( q R)P ( q ), 



(5.6) 



where we used the fact that the function M^(k) defined in equation (|2.2p asymptotically behaves like Mq oc k 2 (k°) in 
the limit of k — > (k — > oo). Here, the quantity cr{R) represents the rms amplitude of the linear density fluctuation 
5q. Then, from (15 .4[) and (15. 5p , we obtain 



b{k; z) ~ b\ 



12 f b2 n<~) g2 ^) 



1/2 



(5.7) 



The expression (I5.7|) contains a term inversely proportional to M^(fc), which leads to the strong scale-dependence at 
k — > 0. Now, we recognize the fact that the scale-dependence of the biasing factor depends not only on R and /nl, 
but also on the ratio of biasing coefficients, b 2 /bi, and redshift, z. 

It is interesting to note that similar kinds of behavior can be found in the halo biasing prescription 0, EH, • 
While the halo biasing describes the clustering properties of virialized objects based on the Press-Schechter theory, 
the origin of the scale-dependence is qualitatively the same as that of the local biasing mentioned above. In this sense, 
the scale-dependent biasing may be a unique character of the local model of primordial non-Gaussianity and it would 
be an important observable of the large-scale structure. 

Note, however, that the scale-dependent biasing factor in the local biasing seems problematic because of the loga- 
rithmic divergence at R — > 0, which originates from the rms amplitude a(R) in the case of a CDM power spectrum. 
Further, the biasing factor b(k; z) can be ill-defined for the negative /nl or b 2 and it eventually becomes complex 
at k — ► 0. McDonald [4(| proposed a renormalization treatment to remove these difficulties by adding a non-local 
counter term to the relation (|5.ip . Although the addition of a counter term might be a possible solution, uniqueness 
of this treatment seems is unclear, and thus the physical reason for the non-local counter term is unclear. In fact, 
as shown in next subsection, the scale-dependence of the biasing factor turns out to be model-dependent, and the 
logarithmic divergence disappears when we consider the equilateral model of primordial non-Gaussianity. 



B. Equilateral model 

In Figure [71 we plot the ratio of the biasing factor, b(k; z) /b G (k; z), in the equilateral model. Here, instead of 
showing the dependence of smoothing radius, we neglect the effect of smoothing (i.e., R — > 0) and just examine the 
effect of scale-dependent non-Gaussianity by changing n to (left), —0.3 (middle) and 0.3 (right). 

In marked contrast to the local model, it turns out that the biasing factors 6(fc; z) in the equilateral model are almost 
constant on large scales, and the linear deterministic biasing, P ga i(fc) = fri-Pmass(fc)) is an excellent approximation to 
the galaxy power spectrum on large scales. These features are irrespective of the values of k, and can be deduced 
from the asymptotic behavior of the function F . In the k — > limit, we have 

1 8 

F(k;z,fa L \R = 0) ~ —f NL D 3 (z)M c (k) 
5 



(5.8) 



Recalling the power-law nature of the primordial spectrum P^(q) oc g™ s ~ 4 , the integral in the above equation converges 
if n s ~ 1 and < 0.3. As a result, the function F scales as F(k; z, /nl|-R = 0) oc M^(k) oc k 2 and the non-Gaussian 
contribution becomes negligibly smaller than the linear biasing term, b 2 P mass (k; z). 

Therefore, as long as the equilateral model of primordial non-Gaussianity is concerned, it is hard to constrain /nl 
through the scale-dependent biasing factor. Although this conclusion comes from the assumption of the local biasing 
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prescription, this would generally hold for other biasing schemes, including the halo biasing prescription. On the 
other hand, a big difference in the scale-dependent biasing factor between the local and equilateral models implies 
that the galaxy biasing is very sensitive to the shape of the primordial bispectrum. An essential reason for having 
a scale-dependent biasing only for the local-type primordial non-Gaussianity is that the amplitude of the primordial 
bispectrum with squeezed configuration, i.e., B^(—k,q,k — q) with \k\ <C \q\, is greater than that of the equilateral 
case (see Eq. (|5.3[) ). The local- type primordial non-Gaussianity has fairly strong mode-correlations between small and 
large Fourier modes, and thereby the biasing, as a small-scale phenomenon, affects the power spectrum on very large 
scales. This generic feature might be very helpful to discriminate between the types of primordial non-Gaussianity. 
The detection of scale-dependent biasing may have interesting implication for the inflationary model, according to the 
single-field consistency relation [l||. The consistency relation states that as long as the inflation dynamics is driven 
by a single scalar field, the primordial bispectrum in squeezed configuration becomes very small, irrespective of the 
type of non-Gaussianity. Thus, single-field inflation generically predicts scale-independent biasing on large scales. 
Therefore, if the presence of scale-dependent biasing on large scales is observationally verified, this strongly disproves 
the consistency relation, and the inflation dynamics may not be simply characterized by the single scalar field. 

VI. DISCUSSION AND CONCLUSIONS 

In this paper, we have comprehensively studied the effect of primordial non-Gaussianity on the power spectrum of 
large-scale structure. Using perturbation theory, we calculate the leading-order non-Gaussian corrections to the matter 
power spectrum through the non-linear mode coupling of the gravitational evolution. In the weakly non-linear regime, 
the signature of primordial non-Gaussianity in the matter power spectrum can be characterized by the primordial 
bispectrum. Adopting two representative models of the bispectrum (local and equilateral), we quantitatively estimate 
the non-Gaussian signals on the matter power spectrum. The primordial non-Gaussianity systematically enhances or 
suppresses the non-linear growth of the power spectrum amplitude at the < 1-2% level. 

We then explore the potential impact on the cosmological parameter estimations from future dedicated surveys 
for BAO measurement. Under the assumption of scale-independent linear biasing, the Fisher-matrix analysis reveals 
that while it is hard to detect the primordial non-Gaussianity to the level of current constraints, the inclusion of 
the non-Gaussian parameter /nl significantly degrades the constraints on the primordial spectral index n s and the 
running of the index a. In this respect, the CMB prior information for n s and a as well as the non-Gaussian parameter 
is very crucial. On the other hand, determination of the distance scale Da{z) is rather insensitive to the presence or 
absence of primordial non-Gaussianity, and thus the characteristic scale of BAOs is a robust standard ruler. 

We have also considered the effects of primordial non-Gaussianity on the galaxy biasing. In the framework of local 
galaxy biasing, in which the number density of galaxies is described by a local function of the mass density field, we 
found that the non-linear mapping of galaxy biasing can modulate the power spectrum and this sensitively depends 
on the shape of non-Gaussianity. In the local model of primordial non-Gaussianity, mode-correlations between large 
scales and small scales are fairly strong, and the non-Gaussian correction induces the strong scale-dependent biasing 
on large-scales, while the scale-independent linear biasing is preserved to a good accuracy in the equilateral model. 
These remarkable properties may be very helpful in discriminating between the types of primordial non-Gaussianity, 
especially in connection with the single-field consistency relation. However, the biasing factor in the local model 
exhibits an apparent divergence, suggesting that the local biasing prescription may be incompatible with the local 
model of primordial non-Gaussianity. Indeed, the local biasing scheme is just a phenomenological parametrization 
and the validity of this prescription itself has not yet been tested enough in the presence of non-Gaussianity. Further 
study using simulations is necessary for quantitative prediction of the non-Gaussian signature on galaxy biasing. 
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